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Abstract 

We study the statistical properties of a cellular automata model of traffic flow with the look- 
ahead potential. The model defines stochastic rules for the movement of cars on a lattice. We 
analyze the underlying statistical assumptions needed for the derivation of the coarse-grained 
model and demonstrate that it is possible to relax some of them to obtain an improved coarse- 
grained ODE model. We also demonstrate that spatial correlations play a crucial role in the 
presence of the look-ahead potential and propose a simple empirical correction to account for 
the spatial dependence between neighboring cells. 
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1 Introduction 



The derivation of coarse-grained descriptions from microscopic dynamics has been an active area 
of research for many decades and a vast amount of literature exists addressing this issue in various 
contexts. In this paper we investigate a particular setup relevant for car traffic modeling. Car 
traffic models can be roughly divided into the following categories (see review papers [SlEllE] and 
references therein): (i) microscopic car-following models that treat cars as particles and postulate 
ordinary differential equations (possibly with delay) for the car velocity; (ii) microscopic discrete 
lattice models where a particular cell configuration with values 1 (car is present) and (car is 
absent) combined with explicit rules for movement between cells are used to represent traffic flow; 
(iii) macroscopic partial differential equations models that are typically conservation laws relating 
the car density and flux; (iv) mesoscopic kinetic models for the velocity distribution, the moments 
of which gives the macroscopic car density and flux. 

While car-following models represent a more realistic setup allowing for very detailed interaction 
rules, lane changing, switching in behavior, etc., lattice models are simpler to implement and are 
more amenable to analytical investigation. Therefore, lattice models have been widely used to 
represent various physical phenomena, including car traffic [HfTTll2"0]. In addition, such cell models 
are closely related to the Ising-type models and cellular automata models and, therefore, a vast 
literature exists addressing various analytical and numerical techniques for models of this type. 

One possible use of vehicular traffic models is to utilize filtering techniques to predict the future 
traffic state. In this context, particle filtering and Kalman filtering for microscopic models has 
been developed (see e.g. [l6l[2U[23] ) • Typically, although very detailed rules of interaction can be 
implemented for microscopic models, Monte-Carlo simulations are necessary to estimate statistical 
properties of car traffic using these models. Such simulations can be very costly computationally 
and using Ensemble Kalman Filtering in conjunction with coarse-grained models is an attractive 
computational alternative [19^124] . 

Considerable effort has been devoted to provide the justification of the macroscopic description 
including derivations of macroscopic models for coarse-grained quantities (e.g. car density) from 
more basic microscopic models. Recently, a novel look-ahead potential was introduced to model 
long-range interactions in prototype lattice model, which were then coarse-grained to obtain macro- 
scopic descriptions (T0Hl2] . In particular, in [20] the look-ahead potential was used to model the 
effect of long-range traffic conditions and a new macroscopic PDE model with non-local interactions 
was formally derived. Extensions to multi-lane traffic have also been developed [HE]. 

In this paper we examine the statistical behavior of the cellular automata (CA) model introduced 
in [20] and the basic underlying assumptions used in the derivation of the subsequent coarse-grained 
model. We demonstrate that while some aspects of the derivation can be improved, assumptions 
about the approximate independence of the neighboring cells does not hold in general. We also 
propose a modification of the macroscopic description based on the numerical evidence for the 
statistical behavior of the microscopic dynamics. Although the new modified macroscopic model 
is empirical, we demonstrate that it can reproduce the results of the stochastic simulations with 
remarkable accuracy. 

The rest of the paper is organized as follows. In section we introduce the CA model, discuss 
the assumptions about its statistical behavior, and following [20], outline the derivation of the 
mesoscopic and macroscopic models. In section [3] we described an improved mesoscopic model that 
is derived by replacing one of the assumptions in [20] with an exact calculation. In section U] we 
present detailed numerical experiments to explore the statistical behavior of the microscopic CA 
model in various parameter regimes, and in section \5\ we discuss the empirical correction to the 
coarse-grained macroscopic PDE that accounts for the spatial correlations at the microscopic level. 
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Conclusions are are presented in section [6j 



2 Cellular Automata Model 

In this section, we summarize the cellular automata model given in [20]. The model describes a 
single class of cars that move in one direction along a single-lane highway. Possible multi-lane and 
multi-class extensions are considered in [HE] • The modeled is defined over a one-dimensional lattice 
C = {1, . . . ,N} of N > 1 evenly spaced cells, and the state of the system is given by a function 
ff:£xM+4{0, 1}^. For t G R + and k G C , 

1 1 , if cell k is occupied at time t ; 
a k (t) = < (1) 
I , if cell k is not occupied at time t . 

Cars are assumed to move from left to right, and only one car is allowed to occupy a cell at a time. 
Periodic boundary conditions are imposed so that o m N+k = 0"fc for any k G C and any integer m. 

Transitions in the state of a are the mechanism for modeling car movement. They obey the 
rules of an exclusion process [Uj: two lattice sites exchange values in each transition and cars may 
not move into occupied cells. In addition cars are only allowed to move one cell to the right. Thus 
the only possible configuration changes are of the form 

{a k (t) = 1, a k+1 (t) = 0} -> {a k {t + At) = 0, a k+1 (t + At) = 1} . (2) 

In the simplest case, the transition rate cq > is a constant. Thus the configuration change in ([2]) 
occurs with probability coAt + o(At) for small At, and in the absence of other forces, cars move 
with velocity vq = hco on average, where h is the cell width. As the car velocity does not depend 
on the cell width, cq must scale with h^ 1 . 

In [20] , the transition rate is given by an Arrhenius-type formula [13] with a look-ahead potential 
J k that depends on values of 07 for I > k. In this case, the transition rates becomes coe~ /3Jk , where 



1 M 

Jk{t) = —^cr k+ i + i{t), (3) 



/3 > is a parameter describing the strength of the look-ahead interactions, and M is the number 
of cars to the right of cell k which affect velocity of the car in cell k. The term e~@ Jk plays the role 
of a slowdown factor when the forward car density is high, i.e, when the road is congested. It is 
also possible to introduce weights into ([3]) so that cars nearby have a more pronounced effect than 
those that are further away. 

The microscopic model described above can be formulated as a continuous-time Markov chain 
and the generator of this process can be computed explicitly. It is defined as 

= lim E*W + (4) 
v YJV ' At^O At 

where ip is any test function and the expectation is taken over all possible transitions between time 
t and t + At. For an arbitrary test function ip, the generator is given by 

= coe- pj "a k [l - a k+1 ] [^ k+l ) - ^(a)] , (5) 
kec 
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where a k,k+l denotes the lattice configuration a after an exchange between the cells k and k + 1: 

(6) 



a 



k,k+l 



ai, l^k,k + l, 
Cfe+i i I = k , 
ya k , I = k + 1 . 

The terms &k an d [1 — (Jk+i] in © ensure that only configuration changes of the form given in (|2|) 
contribute to the sum, i.e., that before a transition can occur, there must be a car in cell k and cell 
k + 1 must be empty. 

A coarse grain model is derived by first computing the action of the generator on the test 
function ip(a) = 07. In this case, the only indices that contribute to the sum in © are k = I — 1 
and k = I. Thus, Ao[ becomes 



Aa. 



1 = c e 



c e 



L o"2-i[l - 07] 
L o"2-i[l - 07] 



l-U 



+ c e-P Jl a l [l-cj l+1 ] 



a 



"I 



c e 



-PJi 



ai[l - a w 



From the definition of the generator, 



d_ 

dt 



(7) 



(8) 



where E is the expectation operator with respect to the probability measure associated with the 
stochastic process. Thus the evolution equation for := E(<Tfc) is given by 



-p fc = E(c t) e-^-i<7 fc _ 1 (l-<7 fc ) 



c e 



(9) 



A mesoscopic model for p = {p\, . . . , p } requires a closure that approximates the right-hand 
side of © by a function of p. In [20], a closure was found under two assumptions: 



AI. The look-ahead interactions are weak so that Ee ^ Jk 



-/3EJ k 



A2. The probability measure on a is approximately a product measure, i.e., f(dk = 8k, ai = Si) = 
f((Jk = dk)f(&i = Si), where / is the probability density. 

Note that Assumption A2 implies, among other things, that and 07 are uncorrelated, i.e. 

E (a k (t)ai(t)) w E(<Tfc(t))E (a^t)) = Pm (10) 

for k ^ I. While this weaker condition is sufficient for closure without the look-ahead, independence 
of higher-order moments is needed for closure with the look-ahead. 

Based on the assumptions above, the resulting equation for the density is 



d 

—p k = c e 
dt 



-pik- 



Vfc-i(l - Pk) ~ c e /3/fc Pfc(l - Pk+i) 



(11) 



where 



1 M 

-Y 

M *H 



Pk+i+l 



and PmN+j = Pi for all integers j. 
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The validity of the two assumptions above depends on the strength of the look-ahead potential. 
Assumption Al implies that (3 is small, in which case the influence of the look-ahead potential Jk is 
weak. Assumption A2 means that (7 kit) and Ok+iit) are independent for all i, including the nearest 
neighbors i = 1,2. This is somewhat unrealistic even without the look-ahead potential, since the 
exclusion principle couples the neighboring cells. Nevertheless, our simulations show that without 
the look-ahead potential the coupling is weak in most situations and the approximate independence 
assumption is quite plausible. On the other hand, in the presence of a strong look-ahead potential 
the neighboring cells become very highly anti-correlated. This can be understood as follows: If the 
length of the look-ahead potential is M, then a car in position k is very likely to be trailed by M 
zeros since the probability to move for the car in the position k — M — 1 is c§e~^l M At <C coAt for 
large /3. Therefore, the probability to move is extremely small, and a car in the position k — M — 1 
is very likely to "wait" for the car in the position k to move. Positive correlations may also occur 
in some situations with and without the look-ahead potential. We discuss such situations in the 
numerical experiments in section HI In summary, assumptions Al and A2 are plausible if the look- 
ahead interactions are weak (/3 <C 1), but are unlikely to hold if the look-ahead interactions are 
stronger (/3 = O(l)). 

A first-order PDE model can be formally derived from the mesoscopic model in the limit N — > 
oo. Let the look ahead distance L and domain length D be fixed; set h = D/N and M = [L/D]N. 
For k £ {1, . . . , N}, let Xk = ik — 0.5)/i denote the center of cell k, and let p N be a smooth function 
of x and t such that p N ixk,t) = Pkit)- Define 

F N (x,t) = e-^ N ^p N ix-h,t)[l-p N {x,t)} . (12) 

where I N is a smooth interpolant of Ik- 

1 M 1 f L 

^)=uE^ + ( i + 1 )'') = T p N ix + y,t)dy + Oih). (13) 
i=i Jo 

Then with a Taylor expansion of F N around Xk, (llip becomes 

§- t P N ix k ,t) = c (F N (x k - 1 ,t) - F N (x k ,t)) = vo-^F N (x k , t) + 0(h) . (14) 

One can then formally pass to the limit — > oo while L and D remain fixed. For each fixed 
x G (0,£>), let k = [x/h] so that Xk — y x as h — y and p — lim^Y— >oo 

p N satisfies 

d d 

-p( x ,t) + —ct>(p)=0, (15) 
where the non-local flux cj) is given by 

<j>{p){x,t) =v e-T^~P^+y^ d yp( x ,t)il-pix,t)). (16) 



3 A New Mesoscopic Model 

In this section we show the expectation of e~@ Jk can be computed exactly so that assumption Al can 
be removed. This leads to a slightly improved mesoscopic model for the car density. (However, as 
in the previous mesoscopic model, assumption A2 is still required to close products.) In particular, 
A2 allows us to approximate 

/ M \ M 

E (e~P J A = E I Yl e~' 3 ' ,Tfc +'+ 1 ) ~ JJe [er^^+A . (17) 

\i=l J i=l 
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where /3' = f3/M. We then use the fact that for 1 < I < N and for any positive integer m, 
{ai) m = oi. Simple algebra gives 



m=0 



00 1 oi\m 



m = 1 



ml 



l + otie-P -1) 



Thus the expectation of e @' ai is 

E( e -^) = 1 + E {at) (e-f - 1) = 1 + pi{e^' - 1) . 
Substituting (|19p into the right hand side of ()1T|) gives 



E e 



n[ 

i=l 



1 + Pk+i+i e 



(18) 



(19) 



(20) 



This expression can be used to derive a new mesoscopic model for the car density pk(t): 



d M 
-^Pk =c pfc_i(l - pk) \\ [l + Pk+i (e _/3 ' - 1 

i=l 



M 

c pfc(l - pk+i) IX [i + Pfc+i+i fe^' 3 ' - 1 

i=l 



(21) 



which can be rewritten as 



Pk =CQPk 



(22) 



One can easily recover the model in (jlip by expanding the exponential and logarithm in ([22 
to leader order in ft . Indeed, simple Taylor expansions give 



log 



1 + Pk+i+l ( e~P' - 1 



'2> 



In particular, for fixed (3 and large M, 

M 



i=l 



l+Pk 



-/3V fc+ i + i + 0(/3 



-/3/ fc + 0{M~ l ) 



(23) 



(24) 



and since M 1 = 0(/i) as iV — >• 00, the new mesoscopic model has the same PDE limit as the old 
one. 



4 Numerical Experiments 

In this section we perform a detailed comparison between the stochastic traffic model and the two 
mesoscopic models (|11|) and (|22p in various parameter regimes. The new mesoscopic model in (|22p 
provides a better front tracking in some regimes. Nevertheless, as discussed earlier, it still suffers 
from the same limitation as the old model, namely that look-ahead interactions should be weak so 
that the measure associate with a is nearly a product measure. When the look-ahead interactions 
are weak {fi is small), both mesoscopic models perform similarly. 
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4.1 The Stochastic Algorithm 

We use the Metropolis algorithm [15] to simulate the stochastic cellular automata model0 Given 
a time-step 5t, we advance the configuration of the lattice during each time-step according to the 
following algorithm. If a car is present at the position k (i.e. a kit) = 1) and there is no car at 
the position k + 1 (i.e. cr^. +1 (t) = 0), then the car at the position k can advance with probability 
coe~" Jk 5t, where J k is the look-ahead potential in (|3]). Simple pseudo-code for a single time step is 
given below. 

Algorithm 1 Metropolis Algorithm (one time step) 
Require: N, M, t, f3, St, and <r fe (i) for k = 1 : N 
for k = 1 : N do 

if a k {t) == 1 AND o-fc+i(t) == then 
J k = Jj Y,i=i a k+i+l(t) 
p = unif(0,l) 
if p< C e-P J *>8t then 
Okit + 5t) = 
<j k+1 it + 5t) = 1 
end if 
end if 
end for 



The configuration array ak(t + St) is stored separately from 0fc(i). Thus, changes in the con- 
figuration in one part of the lattice do not affect the calculation of J k for other indexes in the 
lattice. 

All numerical simulations presented in this paper are ensemble Monte-Carlo simulations and 
expectations plotted in all Figures are ensemble averages. Given the initial conditions (j26|) . we 
generate n = 5000 different realizations of a and estimate the density by the formula 

n 

where the integer p is the realization index. During the course of the simulations, various statistics 
are recorded in order to make the comparisons presented below. 

We perform all experiments using the deterministic "red light" initial condition 

f 1 20 < k < K, 
a J0) = < - - > 26) 

V ' [0 K < k < N, V ; 

with K = 60. This condition allows us to examine the behavior of the leading rarefaction wave. 
Many other configurations are possible including those which lead to shock waves. Following |20j . we 
set the cell size to 22 feet. Other parameters in the simulation are cq = 1/0.23 = 4.3478 (corresponds 
to approximately 65 miles/hour) and N = 700. The boundary conditions are periodic, but N is 
set large enough to ensure that they do not affect the simulations. 



The Metropolis algorithm was cross- validated with a Kinetic Monte-Carlo algorithm [22]. Several comparisons 
were made to ensure statistically accurate results. 
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4.2 Assumption Al: Expectation of the Exponential 

We first check that the new mesoscopic model removes the error in the closure assumption Al, 
i.e. that E(e~^°" fe ) ~ e _ ^ E<Jfe . Two representative cases are given in Figures [1} one with ft = 0.5 
and M = 1 and a second with ft = 3 and M = 5. For reference, we include the formula in (|19D . 
which is used in the new model. As expected, this formula is exact for all values of ft and M. The 
dependence of the results on ft is generic: Assumption Al is nearly valid for values of ft < 1, since in 
this case the powers of ft quickly decay in the expansion of the exponential. However, the accuracy 
of the approximation quickly decreases for larger values of ft. Many other parameter values have 
been tested to confirm this. Roughly speaking, the validity of Assumption 1 appears to be largely 
independent of the size of M. 
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Figure 1: Comparison of E(e /3(7k ) with the closure E(e P° k ) ~ e /3E(7k the closure in (fT9|h Note 
that the solid blue and the red line overlap completely. 



4.3 Front Tracking 

In our next set of experiments, we investigate how well each of the mesoscopic models tracks the 
rarefaction front. We consider the look- ahead distance M = 5 and vary the strength of the look- 
ahead interaction j3. A comparison of the stochastic simulation and both mesoscopic models is 
given in Figure [5J In all cases, the bulk of the traffic proceeds faster in the stochastic simulation 
and the trailing front is less steep, with the new model always outperforming the old one. As 
expected, the difference between the stochastic and mesoscopic models increases as (3 is increased. 
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Indeed, for larger values of /3, there is a considerable discrepancy between the stochastic simulations 
and both mesoscopic models. However, front tracking with the new mesoscopic model f[22l) is still 
considerably better than with the old model (fTT|) . 
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(f) P = 5, M = 5 



Figure 2: Comparison of Simulations of the Stochastic and Mesoscopic Models with the initial 
condition (|26j) . Note that for f3 = 0.5, the numerical results from the mesoscopic models (jlip and 
(|22p overlap in this regime. 
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4.4 Assumption A2: Correlations 

To assess the source of the discrepancy between the microscopic model and the new mesoscopic 
model in (|22|) . we investigate numerically the validity of assumption A2 — that the measure associ- 
ated with a is (nearly) a product measure. To do this, we examine for small i > 1 the size of the 
sample Pearson correlation coefficient [6jH 



where Sk(t) is the standard deviation of the stochastic process at cell k and time t, bars denote 
sample means, and p indexes the samples. 

We compute correlation coefficients for 1 < i < 4 using three sets of parameter values 

1. P = and M = 0, i.e, no look ahead potential (Figure [3]); 

2. (3 = 3 and M = 1 (Figured]); 

3. = 3 and M = 2 (Figure 
The main conclusion are as follows: 

• Without the look ahead, the correlations are relatively weak, except at the trailing front, 
which shows a moderately strong positive correlation (see columns 2-4 of Figure [3]). At 
early time this correlation is stronger for smaller i, but comparable at later times for i > 1. 
These correlations occur because any car that is slower than the "main pack" (maximum of 
the density in Monte-Carlo simulations) will cause cars behind it to also slow down. Since 
it is very likely that more than one car is slower than the main pack, it produces positive 
correlations at the trailing front. 

• When P is nonzero, the correlations rk,k+i depend strongly on the value of i relative to M. 
For i < M, the correlations are large and negative, with the strongest correlation being at 
the trailing front. This is due to the look-ahead potential which tends to maintain a spacing 
of M zeros between cars. For i > M, the correlations are similar in the zero-potential case, 
especially at longer times. 

2 Note that a zero correlation coefficient is a necessary but not sufficient condition to have a product measure. 
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Figure 3: Left scale: correlation coefficients r^+p, 1 < p < 4, for (3 = and M = 0, i.e., no look 
ahead potential. Right scale: mean car density, 
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Figure 4: Left scale: correlation coefficients rk t k+p, 1 < P < 4, for /3 = 3 and M = 1. Right scale: 
mean car density, 
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4.5 Testing the Entire Closure for the New Model 

Next, we test the closure for the right-hand side of the mesoscopic model in ([9]) with the look- ahead 
dynamics, but just one car in the look-ahead potential. We select parameters M = 1 and j3 = 3, 
and using sample statistics, test the closure 



E 



co 



c e 



+1 cr fc _i(l - a k ) - c e 



+2 



(28) 



l + Ecr fc+ i(e-' 3 -i) Ea fe _i(l-Ea fe )-co 1 + Ea fe+2 ( e - /3 - 1) Ea fe (l - E<7 fc+ i). 
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We compare these results to the closure which assumes no look-ahead potential, i.e., we test the 
closure 



E [c e-^ <Tfc + 1 a fe _i(l - a k ) - c e-^ k + 2 a k (l - a k+1 ) j 

(29) 

c E<Tfc_i(l - Ecrfc) - c Eo- fc (l - Ecifc+i). 

Results from these experiments are given in Figures [6] and where we have averaged the 
expressions over five cells to reduce the noisy fluctuations due to the stochastic nature of the 
equations. It is clear in from Figure [6] that the closure in (|28|) is more accurate near the trailing 
front. However, a close inspection of the leading front (see Figure (J7J)) shows that the closure in 
(f29|) does a better job at the leading front. Indeed the closure in (|29|) underestimates the exact 
value; consequently not enough cars move with the leading front. This difference is quite subtle, 
but is observable for times well beyond those give hereJl Moreover, the effect in the density is quite 
noticeable. From Figure ([8]), one may observe the general trend that the closure in (|29p leads to 
better results at the leading front, while the (j28[) performs better near the trailing front. 



The difference becomes lost in the noise around time t = 100. 



15 



Time = 1 



Time = 7 



Time = 14 





I 






RHS 

p=3.0 


) 0.2 0.4 0.6 
Time = 1 





f 


■ 




I RHS 

p=o 


) 0.2 0.4 0.6 
Time = 1 







04 





02 







-0 


02 


-0 


04 


-0 


06 











RHS 

p=3.0 






04 





02 







-0 


02 


-0 


04 


-0 


06 



1 




0.8 




0.6 




0.4 




0.2 














0.2 0.4 0.6 
Time = 7 









l> 

i 


RHS 

p=o 



0.2 0.4 0.6 
Time = 7 




0.04 
0.02 


-0.02 
-0.04 
-0.06 



RHS 

p=3.0 



0.2 0.4 0.6 
Time = 14 



0.04 
0.02 


-0.02 
-0.04 
-0.06 



S 



RHS 

p=o 



0.2 0.4 0.6 
Time = 14 



1 




0.8 




0.6 




0.4 




0.2 















Figure 6: Comparison of closures from (I28j) (top line) and (]29j) (middle line) during the course 
of a stochastic simulation. Data is averaged over five cell intervals to remove noise. The bottom 
line is the density profile, which is provided as a reference for locating the fronts. The look ahead 
potential is M = 1. 
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Figure 7: Comparison of closures from (j28[) (top line) and ([29]) (middle line) during the course 
of a stochastic simulation. Here we plot the differences between the exact closure and the two 
approximate closures. Thus positive values mean and overestimate of the time derivative. Axes 
are magnified to allow inspection of the leading front. Data is averaged over five cell intervals to 
remove noise. The look ahead potential is M = 1. 
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Figure 8: Comparison of the car density in the Monte-Carlo simulations of the stochastic model 
with /3 = 3, 6 M = 1, 3 and simulations of the new mesoscopic model and the mesoscopic model 
without /3 (i.e. (3 = 0). 
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5 Empirical Correction to the Equation for the Density 



As observed in the previous section, the proposed closure in (|28p underestimates the right-hand 
side of (|9|). This is also evident from the plots in Figure [8] comparing the Monte-Carlo simulations 
of the stochastic problem with simulations of the mesoscopic equations. In particular, the bulk of 
cars behind the leading front not propagate fast enough in the mesoscopic model. This means the 
effect of the look-ahead dynamics in the proposed closure (|2Hj) is too strong. 

On the other hand, the right-hand side in ()29p (equivalent to setting f3 = in the right-hand 
side) agrees much closure with the microscopic model. We therefore conjecture that the cross- 
correlations effectively make the look-ahead parameter, (3, in the mesoscopic model smaller. For 
example in Figure El the leading front in the solution of the stochastic model is "in between" the 
two solutions of the mesoscopic model with and without j3. 

Based on these observations, we compensate for the error in the right-hand side of ([9]) by 
introducing a nonlinear look-ahead parameter /3 = (3pf, where d is an empirical fitting parameter. 
Since < p k < 1 the proposed, nonlinear compensation of the look-ahead dynamics leads to the 
effective reduction of the look-ahead potential. The resulting mesoscopic model is 

d M d 

j t Pk = copfe-i(l - Pk) 11 [ l + Pk+i (e _/3 ' p "+' - l) 

i=l 

(30) 

M 

- c oPk (l - p k+1 ) Yl [l + p k+i+l (e-P' p *+i+i - l) . 
i=l 

Results comparing the Monte-Carlo simulations of the stochastic model and the simulations 
of the modified mesoscopic model from (|30p are depicted in Figure [9] for four different parameter 
combinations of f3 and M. As a reference, these results should be compared to those from Figure 
El At this point, d is a tuning parameter. From experiments, we have found that a "good" choice 
of d is sensitive to the value of M, but rather insensitive to the value f3. 
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Figure 9: Comparison of the car density in the Monte-Carlo simulations of the stochastic model 
(blue line) and simulations of the mesoscopic model in (|30|) (magenta). The exponent d in (|30p 
seems to depend strongly on M, but not on j3. 
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5.1 Continuum Limits 

Following the arguments from Sections [2] and [3] it is straightforward to show that the continuum 
limit for the empirical closure takes the form of the conservation law (|15p . but with a new flux 
given by 

=voe-iti~ pd+1 ( x+v > t)dy p{x,t){l-p{x,t)) . (31) 
As in [20] , one may formally expand the integral in the exponential^ To second order, we find 

\ jf p d+1 (* + y. t)dy ~ P d+ \ X ) + —t^H*) + Y&~ pd+l{x) (32) 

so that, via Taylor expansion, 

e -if L p d+ H*+y,t)dy „ e -pp»H*) A _ ^L^_pd + i (x) _ ^&_pi+i {x) \ (33) 

This gives the following local, continuum model 

p t +(e-^ d+1 p(l - p)) x = f (p(l - p)e-^ + \^%)+^ (p(l - p)e-^ d+1 (p^) xx ) x (34) 

which introduces an additional layer of nonlinearity to the models presented in [20]. We note, how- 
ever, that the proceeding derivation is completely formal. In particular, it assumes that solutions 
are smooth and that the error in the Taylor series expansions are small. 



6 Discussion and Conclusions 

In the paper, we have examined in detail the cellular automata (CA) traffic model introduced 
in [20] , including a host of numerical experiments performed on a rarefaction wave evolving from a 
"red-light" initial condition. The work in [20] included mesoscopic (ODE) and continuum (PDE) 
approximations of the cellular automata model. For the mesoscopic model, two assumptions are 
required. The first of these has been removed in the current paper, leading to an improved model 
at the ODE level, but with the same PDE in the continuum limit. Since the look- ahead potential 
conidered here have also been used in other contexts [9|,[T2],[18] , the improved mesoscopic model can 
lead to significant improvements in other areas. 

The second assumption used in [20] to derive the mesoscopic model is based on the independence 
of spatial neighbors in the stochastic process defined by the CA model. We have shown numerically, 
that such an assumption does not hold for strong look-ahead potentials and moreover, that the 
correlations effectively weaken the effect of the look ahead potential at the macroscopic level. We 
conjecture that this effect is density dependent and introduce an ad-hoc macroscopic parameter 
which scales like a power law in the density. The choice of the exponent in the power law was 
determined experimentally and found to depend strongly on the look ahead distance M. However, 
for a given M, the exponent behaves quite well for different f3 and over long time scales. Finally, 
using this new ad-hoc model, we have derived nonlinear variants of the known, local continuum 
models derived from the non-local continuum model in [20]. (See also Section [2]) 

4 This is done for a general potential in [20]; here we consider only the constant potential j3. 
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We see three avenues for future work. First is the need to explore the generality of the numerical 
results presented here. Certainly a complete repeat of the experiments in the current paper should 
be done for a jam — that is, the movement of faster cars at low density into a slower, higher-density 
regime. For a continuum hyperbolic model, this would correspond to a shock. The second issue 
to explore is the relationship between the exponent d and the look ahead distance M. While the 
existence of such a relationship can be inferred by our numerical experiments, a formula would 
prove quite useful in practice. Finally, for engineering purposes like traffic microscopic models are 
not computationally practical. A possible alternative is to use filtering [2], where the macroscopic 
or mesoscopic models are corrected in real-time using measurements. Such an approach has already 
been investigated, for example, in [I61I2T1I23] . 
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